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Abstract 

In nonrelativistic Quantum Chromodynamics (NRQCD), we study the next-to-leading order 
(NLO) QCD radiative correction to the color-octet J ftp inclusive production at B Factories. Com- 
pared with the leading-order (LO) result, the NLO QCD corrections are found to enhance the 
short-distance coefficients in the color-octet J/V' production e+e" —>■ cc(^Pq^^ or ^PQ^^)g by a fac- 
tor of about 1.9. Moreover, the peak at the endpoint in the J/ip energy distribution predicted 
at LO can be smeared by the NLO corrections, but the major color-octet contribution still comes 
from the large energy region of J/tp. By fitting the latest data of a{e^e~' J/ip + Xnon-cc) 
observed by Belle, we find that the values of color-octet matrix elements are much smaller than 
expected earlier by using the naive velocity scaling rules or extracted from fitting experimental data 
with LO calculations. As the most stringent constraint by setting the color-singlet contribution 
to be zero in e+e^ J/ijj + Xnon-cc, we get an upper limit of the color-octet matrix element, 
(0|O'^/'^[i5j^)]|0) +4.0(0|0-^/^pPo(^^]|0)/m2 < (2.0 ± 0.6) x 10-^ GeV^ at NLO in 



a. 



PACS numbers: 13.66.Bc, 12.38.Bx, 14.40.Pq 
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I. INTRODUCTION 



It is widely believed that the heavy quarkonium production and annihilation decay can be 
described by an effective theory, non-relativistic quantum chromodynamics (NRQCD)jl|. In 
the NRQCD factorization approach, the production of a heavy quarkonium is described by 
a series of heavy quark pair states, which are produced at short-distances, and then evolve 
into the heavy quarkonium at long-distances by emitting or absorbing soft gluons. The 
long-distance NRQCD matrix elements are scaled by the relative velocity v between quark 
and antiquark in the quarkonium rest frame. And the heavy quark pair states at short- 
distances should not only have the same quantum numbers as those of the quarkonium, but 
also have other different quantum numbers. In particular, such heavy quark pairs can be in 
a color-octet state. This is the so called color-octet mechanism. 

The color-octet scenario seems to acquire some significant successes in describing heavy 
quarkonium decay and production. But recently, several next-to-leading order (NLO) QCD 
corrections for the inclusive and exclusive heavy quarkonium production in the color-singlet 
piece are found to be large and significantly relieve the conflicts between the color-singlet 
model predictions and experiments. It may imply, though inconclusively, that the color- 
octet contributions in the production processes are not as big as previously expected, and 
the color-octet mechanism should be studied more carefully. Lots of work have been done 
to investigate the color-octet mechanism in NRQCD for heavy quarkonium production. 

The current experimental results on J / 'ib photoproduction cross sections at HERA seem 



to prefer the N 
singlet piece 



O color-octet predication , rather than the description of the NLO color- 



The D 



jPHI data favor the NRQCD color-octet mechanism for 

J/ip production in 77 J/ipX 

The NLO QCD corrections to J/ip and T production at the Tevatron and LHC were 



calculated including the color-singlet piece 



9, 



lOl, 



ll| and color octet piece pJi]. The NLO 



color-singlet contributions are found to significantly enhance the cross sections especially in 
the large pt region. The NLO QCD corrections to polarizations of J/ib and T at the Tevatron 

n n n 

and LHC were also calculated 12, 03, uM ■ The QED contributions to the production of 

n n 

J/ip there and NLO QCD corrections were also calculated [ISj, lla]. The NLO relativistic 



corrections to J/^jj production at the Tevatron and LHC were considered too [l7j. The 



experimental data of polarizations favor the NLO QCD corrections of the color singlet piece 
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rather than the color octet piece. Recent developments and related topics in quarkonium 
production can be found in Refs. Q,Q,0. 

The charmonium production in e~^e~ annihilation at B factories has also provided an 
important test ground for NRQCD and color-octet mechanism. The large discrepancies of 
J/ip production via double cc (including a hidden or an open charm pair) in e^e~ annihilation 



at B factories between LO theoretical predictions 



21 



22 



23|,[2J, |25| and experimental results 



271] once were chal 



NLO QCD 
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eng. 
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:ing 
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issues but now are largely resolved by hig. 
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35| and relativistic 
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37,S8,S9 



ler order corrections: 
4011 corrections, and 



the results show that the color-singlet NLO corrections (both in as and v) may increase the 
cross section of double charmonium production e.g. e+e^ J/'ipVc by an order of magnitude, 
and indicate that the color-singlet contributions are overwhelmingly dominant in most cases 
and there seems no much room for the color-octet contributions in charmonium production 



in e~^e annihi 
be seen in 41 



ation at B factories (discussions on the light-cone and other approaches can 



42 



43 



ii). 



In the J/ip inclusive production e"'"e" ^ J/ip + X atB factories, there are two color-octet 
processes. One is e^e^ qq + J/ip + X{q = u,d,s) studied in Ref . |45l| . where the light 
quark q (or q) emits a hard gluon which turns into a color-octet ^5*1 cc pair fragmenting into 
a J/ip with soft hadrons. But the short- distance coefficient of this color-octet process was 
found to be negligible at ^/s = 10.6 GeV and can only be important at much higher energies 
than ^/s = 10.6 Gevj45|. Therefore, this process can be ignored at B factories. 



The other color-octet process e+e^ ccl^Plf^ or ^PQ'^')g was studied by Braaten and 
Chen|46|]. Based on the LO NRQCD calculation, they predicted that the J/ip production 
is dominated by the region near the upper endpoint in the J/ip energy distribution, and 



3p(8)^ 



150 MeV. But the measured J/ip 
do not exhibit any enhancement 



the width of the peak near the endpoint is of the order of 
spectra in e"'"e~ annihilation by BaBar 47] and Belle {48 ] 
near the endpoint. 

On the experimental side, the total prompt J ftp cross sections in e"*"e~ annihilation were 
measured to be (Xtot = 2.52 ± 0.21 ± 0.21 pb by BaBar j^^], whereas Belle gave a much 



smaller value atot = 1-47 ± 0.10 ± 0.13 pb 



48| . Obviously, the large discrepancy between 



the two measurements should be further clarified. On the theoretical side, for the J/ijj 
inclusive production cross section, the color-singlet contributions including e"^e" —>■ J/ip+gg, 
* J/ip + cc and e"'"e~ — > J /ip + qqgg{q = u, d, s) at LO in as were estimated to be only 
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about 0.4 ~ 0.9 pb which might imply that the color-octet contribution should play 



an important role in the inclusive J/ip production 45|]. However, the ratio of J/ '0 p roduction 
rate through double cc to that of J/ip inclusive production measured by Belle [26] 

J/i) + cc] 



Rr 



a\e ' e 



J/ip + X] 



0.591°;}^ ±0.12, 



(1) 



are much larger than LO NRQCD predictions 



at LO in a,, the ratio is about 0.2 ~ 0.4 



45 



f only including the color-singlet contribution 



49| . And a large color-octet contribution to 



the J / ip inclusive production would enhance the denominator and then further decrease this 
ratio. So, this became a very puzzling issue. Some theoretical studies have been suggested 
in resolving this problem. Fleming, Leibovich and Mehen use the Soft-Collinear Effective 
Theory (SCET) to resum the color-octet contribution jsol. Lin and Zhu use SCET to analyze 



the color-singlet contribution to e"'"e~ J/i^QQ [51]. Leibovich and Liu sum the leading 
and next-to-leading logarithms in the color-singlet contribution to the J/ ip production cross 
section {52 As a new step, Zhang and Chao find the NLO QCD corrections to e+e^ 
J/ip + cc[31] to be large, and increase the cross sections by a factor of about 2 (using the 
same matrix elements as LO), making the ratio R larg er than the LO results. 



Very recently. Belle reported new measurements 



a(e+e" ^ J/ip + X) = (1.17 ± 0.02 ± 0.07)p6, 



(2) 



a[e e 



j/^ + cc) = (0.74 ± 0.08l[J;[J^)p6 



(3) 



cT(e+e- ^ J/i) + Xnon cc-) = (0.43 ± 0.09 ± 0.09)p6. 



(4) 



The J / ip inclusive production cross sectiongiven in Eq. ([2]) is significantly smaller than that 



given previously by BaBar[47| and Belle 



481]. The double charm production cross section 



given in Eq.([3]) also becomes smaller accordingly. The cross section of J / ip + Xnon cc includes 
the color-singlet contribution of e+e^ J/ip + gg and the color-octet one of e+e^ — *■ 
cc(^Po^'' ^-Pq^'*)^ The color-singlet piece has been investigated by including the NLO 
0{as) correction 34 , 35] and O(f^) relativistic correction [40], of which each contributes an 
enhancement factor of 1.2 — 1.3 to the cross section of e+e~ J/ip + gg. As a result, the 
color-singlet contribution has saturated the observed value given in Eq.(jl]), leaving little 
room for the color-octet contribution. 
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In order to further clarify this problem, it is certainly useful to study the color-octet 
process e^e~ —>■ ccl^Plf"^ or ^Plf'^)g itself and to examine the color-octet effect of the next- 
to-leading order (NLO) QCD correction on the J/ip production at B factories. In this paper 
we will focus on the NLO QCD correction to the cross section and J/tp energy distribution 
in e~^e~ — >■ cc paper is organized as follows. In Section II, we will 

calculate the leading order color-octet production cross sections. In Section III , we will 
calculate the NLO virtual and real corrections. In section IV, we will give the numerical 
results and relations to the color-octet matrix elements. A summary will be given in Section 
V. 



II. LEADING ORDER CALCULATION 

In the NRQCD factorization framework, we can write down the cross section for the J/ijj 
inclusive production as 

da{e^e- J/tp + X) = ^d<T(e+e- cc[n] + X){0\O''/^[n]\0) , (5) 

n 

Here da is the inclusive cross section for cc pair in a color and angular momentum state 
labeled by [n] = 2S'+i^W produced in e"^e~ annihilation. And S, L, J is the spin, orbit, 
and total angular momentum quantum numbers of the cc, and i = 1(8) means that cc is in 
color-singlet (octet) state. The short- distance coefficients are calculable in a perturbation 
series in as- The long-distance matrix elements {On'^) are the vacuum matrix elements of 
four-fermion operators in NRQCD [1,]. The long-distance matrix elements are scaled by the 
relative velocity f ^ 1 of the c and c quarks in the cc rest frame. 

At lowest order in v the only term in Eq. ([5]) is the color-singlet contribution , \n \ = ^S^\ 



which is scaled as Olv"^). The coefficient for this contribution starts at 0{al) 45|, |49|, and 



its contribution is away from t 
center-of-mass energy squared 



le upper endpoint Emax = (s + m'^J^^)/{2^/s), where s is the 



451 ]. The color-octet contributions also start at 0{a^), but 
are suppressed by ^ 0.1 relative to the color-singlet contribution, and they are negligible 
in most of the allowed phase-space at leading order in perturbation theory. However, as it 



is pointed out in Ref. 



461], there is an 0{as) contribution to color-octet production near the 

lc(8) 



upper endpoint, for which the Feynman diagrams are shown in Fig. [H Here [n] is ^Sq or 

3(8) 
J ' 



^Pf'\ and contributions of other Fock-states are suppressed by f ^. The three P-wave matrix 
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e+(k2) 

FIG. 1: Two Born diagrams for e~e'^ cc + g. 
elements are not independent at leading order in f ^, and are related by 

(0|O-^/^[=^pf ]|0) = (2J+ l)(0|C-'/^[='Po^'^]|0)(l + 0{v^)) (6) 



In the LO calculation, we refer to e.g. Ref [46[. Momenta for the involved partic 
assigned as e~(fci)e+(fc2) — >■ cc(2pi) + g^k^). In the calculation, we use FeynArts 
generate Feynman diagrams and amplitudes, FeynCalc [55] for the tensor reduction, and 



es are 
to 



CP c 

5^ 



LoopTools 



561] for the numerical evaluation of the infrared (IR)-safe one-loop integrals. 
da is related to the amplitude of created charm quark pair in a color and angular mo- 
mentum state [n], 

A{e^e- ^ccC'^'Lf){2p,)+g{k,)) 

= $^$^$^(si;S2 I SS,){LL,;SS, I JJ,){3j;3k | 8, a) x 

LzSz siS2 jk 

A{e^e' Cj{pi) + Cfc(pi) + gik^)) {L = S), 

e*^{Lz)^A{e+e~ Cj{p, + q) + Ck{pi - q) + gik^)) {L = P), 

q=0 

where (3j;3fc | 8, a) = ^J2T°- , (si; S2 | SS-^, {LLz]SSz \ JJz) are respectively the color- 
SU(3), spin-SU(2), and angular momentum Clebsch-Gordan coefficients for the cc pairs 
projecting on appropriate bound states. ^(e"''e~ — > Cj{pi + q) + c^ipi — q) + gik^)) is the 
scattering ainplitude for the cc production. We introduce the spin projection operators 
PssAva) as|57| 

PssAp, Q) = S2\SSz)v{pi - q; Si)u{pi + q; S2). (8) 

SlS2 

Expanding PsSzipi^ l) i^i terms of the relative momentum g, we get the projection operators 
at leading term and next-to-leading term of g, which will be used in our calculation, as 
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follows 



PisApuO) = ^f{S.Wi+mc). (9) 
Poo(P,0) = -^75(^1+ "^c). 

III. NEXT-TO-LEADING ORDER CORRECTIONS 

To the next-to-leading order in a^, the cross section is 

O = CTBorn + O-yirtual + O-real + 0{a'^al), (lO) 



where 



4 2s 

da^irtual = 4 ^ X] 2 Re{MBo^^MNLo)dPS2 



dareal = \ ^ ^ \Mreal\^ dP S^. (ll) 

The factor 1/4 is the average over spins of initial states. The factor l/2s is the flux factor. ^ 
means summation over the polarizations of initial and final states. Then we need calculate 
M-NLO- Here dPS2(3) means two(three) body phase space. Then the self-energy and triangle 
diagrams are all corresponding to propagators and vertexes of born diagrams. A half of the 
counter term, self energy, and vertex Feynman diagrams for e~{ki)e'^{k2) —>■ cc{2pi)g{k^) 
are shown in FIGS Other diagrams can be obtained by reversing the quark lines. And the 
box Feynman diagrams are shown in FIG El 



At NLO in as, there are ultraviolet (UV), infrared(IR), and Coulomb singularities. We 
choose D = A — 2e dimension and relative velocity v to regularize UV, IR, and Coulomb 
singularities. The self energy and vertex Feynman diagrams contain UV singularities, which 
should be removed by renormalization. The renormalization constants of the QCD gauge 
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CTl 



CT2 



CT3 



SEl 




VTl 



VT2 



VT3 



FIG. 2: Half of the counter term, self energy, and vertex Feynman diagrams for e (A;i)e+(fc2) 
cc{2pi)g{k3). 





BOXl 



B0X2 



B0X3 



FIG. 3: Three box Feynman diagrams for e {ki)e~^ {k2) — > cc{2pi)g{k3). 



coupling constant gs = y/AnOs, the charm-quark mass m and field ^p, and the gluon field 
are defined as 



9's 



m 



(12) 



Here the superscript means bare quantities and the renormalization constants Zi = 1 + SZi 
with i = g,m,2,3. In the NLO calculation, the precision of the quantities SZi should be 
0{as). We choose Z2, Z3, and Z^ in the on-mass-shell (OS) scheme, and Zg in the the 
modified minimal-subtraction (MS) scheme 

a. 



6Z^^ 



6Z. 



OS 

m 



An 
47r 



1 2 , 47ru2 

— + 37b + 3 In — ^ + 4 



-SCf — 



1 



OS 

3 



a 



SZ. 



MS 



^(/3o-2C.) 
2 An [euv 



, 47ru^ 4' 
7ij + In ^ + - 
3 

+ 0(a^ 



+ 0{a 



■ + ln(47r 



C(a 



sli 



(13) 



where /x is the renormalization scale, 7^; is the Euler's constant and /3o = (11/3)Ca — 
(4/3)Tpn/ is the one-loop coefficient of the QCD beta function, and Uf is the number of 
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active quark flavors. There are three massless hght quarks u,d,s so Uf = 3. The charm 
quark c is not included in the running coupling [sl. In this scheme, we do not need to calculate 
the self-energy on external quark and gluon legs. Color factors are given hj Tp = 1/2, Cp = 
4/3, C A = 3 in 5f/(3),. 

Since we are calculating the NLO corrections to the LO cross section, which is already of 
0{as), we have to employ the two-loop formula for Q;s(/i), 

(^sifi) ^ 1 Pi InL 

where L = In (/iVA^^^^), and I3i = {34/3)Ca^ - ^CpTpUf - (20/3)Ca7>«/ is the two-loop 
coefficient of the QCD beta function. 

Exchange of the longitudinal gluon between external charm quark pair in Feynman dia- 
gram BOXl and B0X2 of F1GJ3] should lead to the Coulomb singularities 1/v, where v is 
the relative velocity in the cc rest frame. For the Coulomb- singular color-octet process, we 
find 

da = {Oi'^),oda^'^ f 1 _ !^ + ^ + 0{al) \ . (15) 



The NLO color-octet matrix element {OiJ^) are proportional to ixas/iQv) [l|]. It is just the 
Coulomb-singular term in Eq. ( |T5|) . So that 

da = {Oil^)NLoda^'^ + ^ + ^^"')) " ^^^^ 

Then the contribution of Coulomb singularity has to be factored out and mapped into the 
color-octet matrix element {On'^). 

The Feynman Diagrams in FIGl3]and VT3 in FIG12] contain IR singularity. LoopTools 
can not deal with IR divergent function and five point function, so we need calculate it by 



hand. We can separate the divergence through the way in Ref. 58] . Then all the IR singular 



parts become Co[— /cs,^!, 0, 0, m] and Co[pi, — pi,0,m, m], which are defined as 

J [2tx)^ [q^ - mf/\[[q + piy - ml\[[q + p2Y - 
Cq[pi, — pi,0,m, m] contains soft and Coulomb singularities, and it will appear in Boxl and 
Box2 that are shown in FIG 131 It can be regularized by D = 4 — 2e space-time dimension 
and relative velocity v. 



Co[Pu-PuO,m,m] = ^^^^^{^-^j 1(1 + e) 



2 



- + — -2 + 0(e) 
e V 



(18) 



REALl 




REAL4 



cc 



REAL2 




REALS 



cc 



REALS 




REAL6 



FIG. 4: Six of the twelve real correction diagrams. 



where v = y— (pic — Pic)'^/'m- In the meson rest frame, we have v = \pic — picl/m, which 
is just the relative velocity v between c and c. Co[— fcs, pi, 0, 0, m] has soft and collinear 
singularities, and it will appear in VT3 in FIGJ2] and Box3 in FIGjS] 



C{-k3, -fcs - p, 0, 0, m) 



26 



-e(7£;~ln47r) 



167r2 



- - - In 

e'^ e 



+2ln' 



+2Li2 



Apr ks 

2pi- k^ + 



-2py k-i 



(19) 



There are twelve diagrams for real corrections, six diagrams for gg process, two for qq 
process, and two for ghost process. Half of them are shown in FigHl The other six diagrams 
can be gotten through reversing the charm quark lines. The calculation of the real corrections 
is similar to the leading order calculation, but there should appear IR singularities. \5^ . 



IV. NUMERICAL RESULT AND COLOR-OCTET MATRIX ELEMENTS 



We now turn to numerical calculations for the cross sections. Taking rrij/^ = 2m, m = 
1.55 GeV, A^=388MeV, then a^f/i) =0.245 for fi = 2m. The cross section at LO in is 

iVlo 



a(e + e" 



J/i^ + X) 



11 







lc(8)l 



GeV3 

while the cross section at NLO in a, becomes 



3p(8)i 



GeV5 



pb; (20) 
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1q(8)i 



GeV3 
10 



35 



3p(8)l 



GeV5 



pb. (21) 



= 10.6GeV m = 1.55GeV 
A3 = 0.388GeV 



o . 



O _ 



NLO 



LO 



_j I I I I I I I I I I I I I I I I I 

1.5 2.5 3.5 4.5 ^(^'eF) 

FIG. 5: The short-distance coefficient cj(^S'q^^) in e^e" cc(^sl^^)g as functions of the renormal- 
ization scale fi. 

S_ -I >c > 

s = lO.eCeV m = 1.55GeV 
A3 = 0.388GeV 



NLO 




ti{GeV) 



FIG. 6: The short-distance coefficient (j(^Pq^^) in e'^e —>■ cc('^PQ°')g as functions of the renormal- 
ization scale /i. 

The NLO short- distance coefficients are larger than the LO coefficients by a factor of about 



^3p(8)^ 



1.9. Our LO result is consistent with that in Ref . [46l| . Dependence of the short-distance 
coefficient o"(^S'q^^) on the renormalization scale is shown in FIGJSl and o"('^Pq^^) is shown 
in FIGJi 

If we choose an energy cut Ecut for the J/ip in the endpoint region, we can get the cross 



section cr\Ej/^>EcuT- The result for cc{^S^^) is shown in FIG. Hand the result for cc{^Pj 



3p(8)^ 
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FIG. 7: The short-distance coefficient d'{^ S(^'')\ej^^>Ecut ™ 
cc energy cut EcuT- 



ccQ-S^^)g as functions of the 



shown in FIG. [HI We see that at LO the short- distance coefficients only contribute when 
the J /if) energy is near the endpoint, and the NLO QCD correction can, to some extent, 
smear the J/ip energy distribution near the endpoint. Nevertheless, at NLO the most color- 
octet contribution still comes from the large energy region, say E > 5 GeV. The evolution 
from the short- distance color-octet cc to the final state J/ should affect the distribution of 
J/ip energy via emitting or absorbing soft gluons. If we ignore this effect, the distribution 
of J/ip energy in e"^e~ center of momentum frame is a delta function 6{E — E^ax) and 
Emax = (s + mjy^)/(2-yi) at LO. Braaten and Chen analyzed this evolution effect, and it 
might broaden the energy distribution in the order of 150 Mevj4^. But it does not change 
the fact that most of the color-octet contributions concentrate on the large energy J/ijj 
region. Experimentally, the observed cross sections do not exhibit any enhancement near 



the endpoint 



The color-octet matrix element 



(22) 



was previously extracted from the Tevatron data for J/ip production with the LO calcu- 
lations, which is listed in TableHl If we use the minimum value of the matrix element, 
s = 4.54 X lO^^GeV'^, we can get the cross section at LO in from Eg. (1201) 



cj(e+ + e" ^ J/tp + X) = (0.50 ~ 0.55)pb, 



(23) 
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^ = 10.6GeV m = 1.55GeV 
A3 = 0.388GeV ^ = 2m 




CUT 



FIG. 8: The short-distance coefficient (t{^Pq )\ej^^>Ecut ™ ~^ cc(^Pq )g as functions of 

the cc energy cut Equt- 

TABLE I: Color-octet matrix elements Mk = {0\O-^/'^[^S^^^]\0) + k {0\O'^/'''[^P^^^]\0) /ml 





k 


Mfc(10-2GeV3) 


Cho and Leibovich [60, ^Ij 


3 


6.5 ± 1.5 


Braaten, et al. MRST98LO[62j 


3.4 


8.7 ±0.9 


Braaten,ei al. CTEQ5L[62] 


3.4 


6.6 ±0.7 


Kramer [63] 


3.5 


4.54 ± 1.11 



and the cross section at NLO in as from Eq.(l2T]) 

cj(e+ + e' J/i) + X) = (0.93 ~ 1.08)pb. 



(24) 



The calculated cross section in Eg. ( 12^ for this color-octet process, which belongs to the 
J/ip + Xnon cc production, is much larger than the latest observed value given in Eq.(j4]): 
a[J/i, + Xnon cc] = 0.43 ± 0.13 pb [53]. 

In fact, for the e~^e~' ^ J/ip + Xnon cc production, the color-singlet process e+e^ — *• 
J/ip + gg has been found to make a dominant contribution to the cross section: a{e~^e~ 



0.4 ~ 0.7 pb at NLO in and v'^ 



34 



35, 



401 ] , thus leaves little room to the color- 



octet contributions. This gives a very stringent constraint on the color-octet contribution, 
and may imply that the values of color-octet matrix elements are much smaller than expected 
earlier by using the naive velocity scaling rules or extracted from fitting experimental data 
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with the leading-order calculations. 

Even if we disregard the dominant contribution of e^e'^ J/ + gg to e^e~ — *• 
[J/'ip + Xnon cc] by setting the color-singlet contribution (j(e+e^ J/ip + gg) to be zero, 
and combining Eq.(j3]) with Eq. fl^Tl) . we can get an upper bound of the color-octet matrix 
element: 

M^^^ < (2.0 ± 0.6) X 10~2 GeV=^ (25) 

All the values of the color-octet matrix elements listed in TabU are larger than this upper 
bound by at least a factor of 2 . 

V. SUMMARY 

In summary, we find that the NLO QCD radiative corrections can enhance the short dis- 
tance coefficient of color-octet J /ip production at B factories via e"'"e~ — cc 
with a K factor (the ratio of cross sections of NLO to LO) of about 1.9. The NLO QCD 
correction smears the J /if) energy distribution near the endpoint. But the most color-octet 
contribution is still from the large energy region. The values of color-octet matrix elements 
are much smaller than expected earlier by using the naive velocity scaling rules or extracted 
from fitting the experimental data with the leading-order calculations. If we ignore the 
dominant color-singlet contribution by setting the color-singlet contribution to be zero, and 
use the color-octet contribution to saturate the latest observed production cross section 
a{e^e~ J/ip + Xnon cc), we get the most stringent upper bound for the color-octet matrix 
element: (0|O-^/^[^Sj^^] |0) + 4.0 {0\O^/'^fPl^^^]\0)/ml < (2.0 ± 0.6) x 10"^ GeV^' at NLO in 
as- 
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